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Abstract 

We consider the problem of a scalar field, non-minimally coupled to gravity through a 
—^(p'^R term, in the presence of a Brane. Exact solutions, for a wide range of values of 
CN ' the coupling parameter ^, for both ^-dependent and (^-independent Brane tension, are 

^2i ' derived and their behaviour is studied. In the case of a Randall-Sundrum geometry, a 

t:j- . class of the resulting scalar field solutions exhibits a folded-kink profile. We go beyond 

^ ! the Randall-Sundrum geometry studying general warp factor solutions in the presence of 

Q I a kink scalar. Analytic and numerical results are provided for the case of a Brane or for 

smooth geometries, where the scalar field acts as a thick Brane. It is shown that finite 
geometries with warp factors that asymptotically decrease exponentially are realizable 
^. for a wide range of parameter values. We also study graviton localization in our setup 

r-| ! and find that the localizing potential for gravitons with the characteristic volcano-like 

profile develops a local maximum located at the origin for high values of the coupling ^. 

X 



1 Introduction 

The idea of realizing our Universe as a defect in a higher dimensional spacetime, although 
not newpP, has received a lot of attention in recent years in the framework of String 
Theory where D-Branes^, i.e. membranes on which the fundamental string fields satisfy 
Dirichlet boundary conditions, play a significant role. In the framework of String/M- 
theory[.3j or the AdS/CFT correspondence [4j, Brane models[H][B'^jH] have revealed new 
possibilities for the resolution of the hierarchy problem of particle physics as well as 
for the relation of gravity to the rest of fundamental interactions. In D-Brane models, 
Standard Model fields are trapped on the Brane, while gravitons propagate in the full 
higher dimensional space (Bulk). In an interesting case of a Brane Model with an infinite 
extra dimension, gravitons are localized on the Brane due to the curvature of the extra 
dimension^^^. A solution to the Einsten's equations of motion with a flat metric on the 
Brane and AdS^ geometry in the Bulk exists, provided the positive Brane-tension is 
finely tuned versus a negative Bulk cosmo logical constant. 

Although the Standard Model fields are assumed to be localized on the Brane, grav- 
ity is not necessarily the only field propagating in the Bulk. A number of Brane models 
with Bulk scalar fields have been constructed ^O]^!], either from a theoretical or phe- 
nomenological viewpoint [12j. Actually, the Brane itself could be a defect substantiated 
by a Bulk scalar field configuration (a "kink")^^. The presence of a Bulk scalar field 
opens the possibility of a direct coupling of this field to the curvature scalar. A specific 
form of this coupling corresponds to the gravitational term appearing in the so-called 
tensor-scalar theory of gravity^!]. A Bulk scalar field non-minimally coupled through a 
coupling of the form 0^ R has also been considered in the Randall-Sundrum framework 
and numerical solutions have been discussed 1151. 



In the present article we consider a 3-Brane embedded in 5D space endowed with a 
Bulk scalar field 0, non-minimally coupled to gravity through a — ^ 0^ R term. We inves- 
tigate analytically the existence of solutions to the coupled system of equations of motion 
for the metric and the scalar field in the framework of a metric ansatz Diag (^C^^^^^f^u, l) • 
In the case of the Randall-Sundrum form of the metric we derive analytically a complete 
set of exact solutions for a range of values of the non- minimal coupling strength ^, cor- 
responding to specific choices of the scalar potential. Scalar fields, with or without non- 
minimal coupling, are often introduced against a given Randall-Sundrum background 
under the assumption that their effect on the background geometry will not be im- 
portant. We do find exact non-singular scalar field solutions compatible with an exact 
Randall-Sundrum background, taking into acount the full back-reaction of the field. 

We show the existence of a class of solutions for a general warp function with an 
asymptotic Randall-Sundrum AdS^ behaviour. In all these considerations we allow for 
a field-dependent Brane-tension. Furthermore, we discuss the existence of smooth AdS^ 
solutions for which the role of the Brane is played by a "kink" configuration of the Bulk 
scalar field itself. Both numerical and an approximate analytic treatment of the problem 
is provided. In particular, we calculate the warp factors for smooth geometries in the 
presence of the kink for different boundary values at the origin and obtain various solu- 
tions. Although we concentrate on Z2 symmetric solutions, smooth asymmetric solutions 



are also possible. Through an analytical investigation we verify that for a wide range of 
values of the parameters, we can get warp factors that decrease exponentially and thus 
provide us with finite geometries. We also find analytical solutions for certain special 
values of the parameters and different ranges of ^. In the final section of this paper we 
study graviton localization in our setup and check the form of the localizing potential 
for gravitons with the characteristic volcano-like profile. We find that for ^ higher than 
some specific value the potential develops a local maximum at the origin, which gradually 
increases as we move towards higher values of the coupling parameter. 



2 The Framework 

Consider a general 5D theory of a real scalar field coupled to gravity. Allowing only for 
terms linear in the Ricci scalar, we may write the general Action as 



S = J d'xV^ lfi<P)R -^iV<pf 



Vi<P)-C^\, (1) 



where /(</>) is, for the moment, a general smooth positive-definite function of the scalar 
field (f). Gmn is the five-dimensional metric, not to be confused with the Einstein tensor. 
In the case of a constant /, we have the Einstein Action. The last term corresponds to 
0- independent matter. Note that, the above Action can always be transformed through a 
conformal transformation Gmn — *■ GMNfi'P)/'^^^ iiito an Action where the Ricci scalar 
enters in the Einstein fashion as (2M^) R. Nevertheless, a 0-dependence will arise in the 
matter term giving a theory different than the one we would get in the absence of /(0). 

The equations of motion resulting from (0) are 
/(0) (Rmn - Igmnr) - VmV;v/(0) + GMivVV(0) = ^TifN + ^^S > (2) 

V^'/>-^ + ^7^ = 0, (3) 



with 



Ti% = Vm^P'^n^ - Gmn ( ^(V^)^ + V(0) ) (4) 



the energy-momentum tensor of the scalar field and TJJ^ the energy-momentum tensor 
of (other) matter. 

At this point we shall restrict the metric Gmn introducing the warped ansatz 

3^(^)r/M- 
Gmn = 5 (5) 



where x^ = {x^, x^) = (x^, y) and ry^,^ is the AD Minkowski metric with signature 
(—1, 1, 1, 1). We can always choose ^(0) = 0. 



The presence of a Brane introduces an extra term 



d^xV^a{(l))6{y) = - / d^xa{(f))6{y) , (6) 



where the Brane tension cr(0) is, in general, 0-dependent. Introducing this term in the 
Action, modifies Tmn in Enstein's equations as 

STmn = -G^Jl^d"^ o-(0) 6{y) = -Vf^u^M^N ^(0) ^iv) ■ 

In what follows we shall ignore the presence of (extra) matter beyond the Bulk scalar 
field. Substituting this metric ansatz into the equations of motion and assuming that the 
scalar field is just a function of the fifth coordinate, i.e. = (f){y), we obtain 

^/(if + 2i/ = 1(0)2 -i\/, (7) 

If A + IfiAr + ^-fA + f = -i(0)2 -^V- ia(0) S{y) , (8) 

+ 2i0 - ^ + I (-4i - 5(i)2) - p{y) = . (9) 

The dot signifies differentiation with respect to the fifth coordinate y. 

The Junction Relations at the point y = where the Brane is located are 

A0(O) = 0(+O) - 0(-O) = j^s'^Jy . (10) 



f + W) 



Ai(o) = i(+o) - i(-o) = - ( ^r-w7^ ) • (11) 



The prime signifies differentiation with respect to 0. By a and a' we indicate the values 
at 0(0). 

Among the above three equations of motion in the Bulk only two are independent. 
They can be written as 

\/(0) = -3/i2 - y-A - ^-Af - f , (12) 

l^' = -lfA+^-Af-f. (13) 

Let us now restrict the coupling function /(0) to be a function quadratic^ in 0. 
Introducing a dimensionless parameter ^ and normalizing it appropriately, we may write 

/(0) = 2M3-|0^ (14) 



^Even for a general coupling function, wc may consider an expansion in even powers of the field 
f{(j)) w /(O) + hf'{0)4P + ■ ■ ■ and retain the lowest non-trivial term. Such an expansion would be valid 
for small field values {(f) « (|2/(0)//'(0)|)i/2). 



The scale M is related to the 5D Newton's constant G as 2M^ = (167rG')"^ With this 
choice, we have 

and the equations of motion in the Bulk become 



2^ __(2M^-l<^^jA-^A00 + ^00 + ^0^ (16) 

The Junction Relations take the form 

^^^ 2M3-|(l-fe)02(o) ' ^''> 

AA = - ^ , , ^ ^ , , . 18 

Again, a and a' are the corresponding values at = 0(0). Note the simplification of the 
denominator at the D = 5 conformal value^ C,c = 3/16. 



3 Randall- Sundrum Metric 

In this section we shall make the definite choice of the warp function A{y) to be the 
standard Randall- Sundrum warp function A{y) = —K,\y\ and impose Z2 symmetry on 
the scalar field {(f){—y) = 4>{y), 0(+O) = — 0(— 0)). Substituting, we obtain the equation 
for (f){y) in the y > Bulk 



The values on the Brane will have to obbey (1 - 2^0^(0) = ^0(0) U0(O) + 20(0) j 0(0). 

Thus, the boundary value 0(0) = is possible only with ^ = 1/2. We proceed distin- 
guishing the two cases (0(0) = and 0(0) 7^ 0). 

3.1 Special case with 0(0) = 

In this case, possible only for ^ = 1/2, we have the solution 

0(^) = M±^ (1 _ e-H/2) . (20) 

K 

The Junction Relations give 

a = 12kM^ , a' = 20(+O) . (21) 



^The conformal value in D dimensions is ^c — 4(j)_A ■ See ref.|16|. 



Note that this solution is possible only with field-dependent Brane-tension. Note also 
that the first is the standard Randall- Sundrum relation. This special solution has the 
shape of a folded kink is plotted in Figure 1. Beyond a small region near the Brane it 
reaches a constant value 0(±oo) = 2k^^(J){+0) (Figure 1). 
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Figure 1: Solutions for a field-independent Brane-tension and different values 

The required positivity of the coupling function 

/(0) = 2M' - ^-<p\y) = 2M3 - 2 (0(O)/«:)' (l - e-^l^l/^)' > (22) 

imposes the boundary value constraint 

(0(0))' < K^M^ =^ {a'f < ^-na. (23) 

o 

The scalar potential corresponding to this solution can be obtained to be 

V{(j)) = -2Ka--a'^ + -{K(f) + a'y. (24) 



In the expression above, we have made use of the Junction Relations. If we were to 
start with a general quadratic potential V^c/)) = A + Cicf) + 6*20^5 the solution ^lU^ and 
A{y) = -f^\y\ is possible for A = -2Ka + a'^/8, Ci = na' /2 and C2 = k^/A. 

3.2 General case 

In the general case 0(0) 7^ 0, the equation of motion can be written as 



1. ^0 ^ </> 

2^ ^^0 2 ^rh 



(25) 



and leads to the solution 



0(y) = 0(0) [1 + a(4 - r') (1 - e-'^l^l/') ] '^-' 



(26) 



where 






The Junction Relations take the form 

^^^.r,. {^Ml^Mi^)y±l^m 2:: \--1m-' (28) 

^ ^ 2M3 - I (1 - f^) 02(0) ' 2M3 - I (1 - f ^) 02(0) ■ ^ ^ 

These two constraints can be rewritten as 

cr' = 2K0(O)(a-4O, a = 6fi; ( 2M^ - ^^^(O) ) + 4^/€a02(O) . (29) 



In order to study whether the positivity of the couphng function and the requirement of 
a positive tension Brane {a > 0) can be simultaneously satisfied, we consider the four 
possible sign choices of the non-minimal coupling strength parameter ^ and the boundary 
values parameter a. 

They are satisfied^ 

1) If ^ > 0, a > 0, always. 

2) If ^ < 0, a < 0, always. 

3) If ^ < 0, a > 0, only if 

3 3M^ 



4) If ^ > 0, a < 0, only if 



3 3M3 



3.2.1 Field-independent Brane-tension (a' = 0) 

To simplify our analysis we may consider separately the case of field-idependent Brane 
tension {a' = 0). In this case the Junction Relations simplify to 

a 

[2M3-|(l-fe)0^(O)] 

and 

a = 4e . (33) 

Notice that the positivity of the Brane-tension is always satisfied, since 2M^ —C,(f)'^{0)/2 + 
8(^^0^(O)/3 is positive if the coupling function is positive (2M^ — ^cf^ /2 > 0). 

The relation a = 4^ simplifies the solution ()26p to 

<P{y) = m [1 + 4(4^ - 1) (1 - e-'^l^l/^) ] ^ . (34) 

^We choose 0(0) > 0. 



For C, > C,c = 3/16, the quantity in brackets 16{C,—C,c) — 16(,^ — ,^c — l/4)e '^'^'^^ is positive. 

Note that for the special value ^ = 1/2, if 0(?/) is the solution, so is (j){y) + const.. 
Also, in the special case ^ = 1/4, the solution takes the form 



0(1/) = 0(0) e 



2(l-e-«l!'l/2) 



(35) 



This is shown in Figure 1. 

The requirement of the positivity of the coupling function, in the allowed range 
^c < ^ corresponds to the inequality 2M^ > §0^(0) [1 +4(4^ - l)]^?^. For the special 
value ^ = 1/4, this corresponds to 2M^ > y0^(O). 

For the conformal value C,c = 3/16, the solution reduces to an increasing exponential 



My) = 0(O)ei'^l^l . 



(36) 



For values of the coupling parameter in the range < ^ < C,c the quantity in 
brackets vanishes at yo = ±Mn [1 + (^c — 0^^]) while the exponent is negative, i.e. 



_1 1 < 0. Thus, in this range the solution is singular. 



For ^ negative, the solution 



0(y)=0(O)[4(l + 4|e|)e-''l^l/2-(3 + 16|e|)] 



2|g| 

4iei+i 



(37) 



is characterized by an exponent between and 1, while the expression in brackets vanishes 
at yo = ±- In [1 + (3 + 16|^|)~^]. Note that 0(|/o) = — oo. As we shall promptly see, these 
"solutions" are not acceptable since the scalar potential, possesing negative powers of 
the scalar field, is singular. 

The scalar potential corresponding to the solutions found can be immediately ob- 
tained from equation ()15|). In order to do that it is usefull to obtain the derivatives of 
the solution. They are 



and 



0(l/) = 0(l/) 



0(y) = -0(l/) 





\ 
v 


r 1-4? 
(1+4(4^-1))^ -1 

V0(o)/ 




'm 


-1) f(t)(y)\ 2« 
2^ V0(O)/ 



+ 1 



(3^ 



(39) 



Substituting the above into (fTK|l we obtain 
V(0) 
where 



-6k^M^ + (j)'^ \Ci + C2 



0(0) 



i-4e 
2e 



+ C. 



0(0) 



l-4g 



^-?1^ + V 



e 



+ 4 



4^-1 



(40) 



(41) 



C2 = -4e«:2 (1 + 4(4^-1)) 

positive in the range < < 
potential includes logarithmic terms. It is 



a 



All powers are positive in the range < ^ < |. In the special case ^ = 1/4, the scalar 



^ 



C2 



^(0) 



-2 /I #3 



Qk^M^ + ^K 



2 J 2 



-i + ^ln(0/0(O)) + i(ln(0/0(O)) 



For the special value S, = 1/2, the scalar potential has the quadratic form 



V{<l)) = -6k^M' + ^K^ (30^ - 



+ 250^(0)) 



(43) 



(44) 



For the limiting conformal value ^c = 3/16, all the above coefficients vanish and we 
obtain a constant potential V = —Gk^M"^. For negative values the appearing powers 0^5 

-—2 

and 0« are negative and, since the solution (j){y) vanishes at a finite point, the potential 
is singular. 

3.2.2 General case with field-dependent Brane-tension (a' ^ 0) 

In order to investigate the behaviour of the scalar field solution ()26|) . we first consider 
the case a > 0. In this case we have a solution increasing near the origin, since na = 
0(+O)/0(O) > 0. The quantity in brackets is positive, provided a(4 — ^~^) > —(1 — 
g-K|y|/2^-i_ xhe lower limit of the right hand side is —1, which corresponds to the range 



^ > jT^-rr. As examples, consider the cases ^ 



4+a 



a 



1 and ^ 



a 



1/8. The first one 



corresponds to a positive exponent 2/3, while the second corresponds to the exponent 
— 1/2. They are both shown in Figure 2. 
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Figure 2: Solution profiles for a field-dependent Brane-tension. 

For values of ^ below this bound there is a point for which the expression in brackets 
vanishes and, since, the exponent is negative, there is a singularity. This point is yo = 
— -^ In [1 — ^/a(l — 4^)]. For the special value ^ = -^^-^ the singularity is pushed to 
infinity and we obtain a purely exponential form for the solution, namely 



0(y) = 0(O)e 



Ka\y\ 



(45) 



Before we move to consider negative values, let us mention again the special value 
^ = 1/4 which corresponds to 



m = 0(0) e^"(^-~'='^'/^) 



(46) 



As we have remarked earlier, for ^ > and a > 0, the positivity of the Brane- 
tension (a > 0) is always true. However, the requirement of a positive coupling 
function introduces a constraint on the parameters. It is sufficient to have 2M^ > 

|0^(O) [1 + a{4 — ^^^)]'^^''^ ■ For the special value ^ = 1/4, this constraint has the form 
2M3 > ^02(0). 

For negative values ^ < (and still a > 0) the solution takes the form 



</.(?/)= 0(0) [l + a(4+|er^)(l-e 



2iei 
-k|j/|/2\1 WFT 



(47) 



Note that although the scalar potential has negative powers (^^^/^'^l and ^^^^-'^/l^l), 
there is no singularity, since the scalar field does not vanish anywhere. Note also that for 
^ < and a > the positivity of the Brane-tension introduces a constraint a < |+ ^^ 



4 ' |?|02(O)- 



Let's move now to consider the case a < 0. Writing the solution as 
0(y) = 0(0) [ 1 - |a|(4 - r')(l - e~-\y\/') ] ^ , 



(4J 



we see that for ^ > 1/4, the exponent is positive. For \a\ < | the quantity in brackets 
stays positive. However, for |q;| > i it is necessary to limit the range of ^ to .^ < ^i . 



For the critical value ^ = ^J"_^ , the solution becomes a decreasing exponential, namely 

0(l/) = 0(O)e-'^HIy| (49) 

As examples of the solution in the above range, let's consider a = —1/8, ^ = 1 and 
a = —1/2, C, = 1/3, shown in Figure 3 




Figure 3: Field-dependent Brane-tension solutions with negative a. 



For \a\ > 1/4 and C, outside of the above range, i.e. C, > 



we obtain solutions 



'R 



that vanish at a finite distance from the Brane, namely \yo\ = — - In [1 — ^/|q;|(4^ — 1)]. 

9 



This, again, corresponds to a singular scalar potential due to the negative power ^^^^^ ^ 
that appears in it. 

For values below 1/4 but positive, the solution takes the form 



0(l/)=0(O)[l + |a|(r'-4)(l 



2« 



-t^\y\/'^\ 1 i-4« 



(50) 



and gives a smooth decreasing profile, just as seen above. 

Before we move to consider negative values of ^, let's consider the special case ^ 
1/4. In this case, we have 



(/.(!/) = 0(0)6 



-2|a|(l-e-'^l!'t/2) 



(51) 



shown for a = — 1 in Figure 4 




Figure 4: Profile of the solution for the special value i = \- 



The positivity of the Brane-tension, for .^ > and a < 0, as was found earlier, 
introduces the constraint 

On the other hand, the positivity of the coupling function, since 0(t/) is a decreasing 
function, is covered by 2M^ > |0^(O). 

For negative values .^ < 0, the solution can be written as 



0(1/) =0(0)= [l-|a|(|er^ + 4)(l 



2|e 
= -«|j/|/2\ 1 4151+1 



] 



(53) 



It is easy to see that for \a\ > 1/4, the scalar field vanishes at a finite distance from 
the Brane, namely \yo\ = — | In [1 — |^|/|q;|(4|,^| + 1)]. This, again, amounts to a singular 
scalar potential due to the negative powers 0~^/l^l and 0~(^/l^l+^). 

The scalar potential has exactly the same form as in the o"' = case {a = 4^), the 
only difference being a slight change in the coefficients C2 and C3 which become 



C, = -A^n^[l + ^{A^-1)](-^ 



C, 



CI 



32^^' 






(54) 



10 



4 Beyond Randall- Sundrum 

Let us consider again our original set of the two independent equations (fT^ . (fT!?jl for 
the specific choice of couphng function /(0) = 2M^ — ^0^/2. If we do not impose any 
restriction on the scalar potential function V^(0), we can consider the first equation as an 
equation that determines the scalar potential in terms of the functions 0(|/) and A{y). 
Concentrating on the second equation, we can view it as an equation for the warp factor, 
giving a different A{y) for every different choice of 0(?/) configuration. Motivated by the 
form of the solutions found in the Randall- Sundrum case, we may start by introducing 
a scalar field configuration in the form of a folded kink with Z2 symmetry'^ 



<P{y) = 00 tanh(a|?/|) 



(55) 



where 0o = a ^</)(+0). For the positivity of the coupling function it would be sufficient 
to require 0^(0) < AM^a'^/^. 

Substituting (|^ into the Junction Relations, we obtain, since 0(0) = 0, 



i(+0) 



a 



12M3 ' 



(+0) 



a 
~2 



(56) 



Note that the positivity of the Brane-tension requires A(+0) < 

In the y > bulk, we have = a0o( 1 — 0^/0o), = 
these into the equation of motion, we can write it in the form 



2a20( 1 - 0V0o)- Substituting 



2M 



3 yA x'(0) + |0o0x(0) + 



'+a<pi[^-^ 



1 



721=0' 



(57) 



2' \ 2' J ''' 2 

where X = A{y) and A = X = 0X'(0). This differential equation can be integrated to 
give 



X(0) = Coil- 



k2 \ 1/3 



AM^J 
where Ci and C2 are given by 



+0 <( Ci ( 1 - 4^] 2Fi{l/2, 1/3, 3/2, e0V4M3) + C2 



4M3 



Ci = -- 



r' - 6 + 



12M3 



(2e-i: 



Co = — 



^"-^+4S^^'^-') 



(58) 



(59) 



and the integration constant Co = ^(+0) = — cr/12M^ should be negative in order to 
have a positive Brane-tension. 

The metric warp factor will be given by e"^'-^-' = exp [J (i?/X(0)]. Near the Brane, 
i.e. for y ^ or — i> 0, we have 



^My) ~ gxp 



a 



12M3 



y + 0{y' 



(60) 



*A(+0) = -A(-O), 0(+O) = -0(-O). 
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which is a pure Randall-Sundrum behaviour. 

In the asymptotic region { y ^ oo =^ ~ 00 = a~^0(O) ), we have 



.My) 



exp 



(a0o 



,-1 



;i-0V0g) 



X(0) 



1 



with 7t = — X(0o) and 



■x( 



6o)X'(0o) + ---) ^e-^^^(i/), 



9 



(y) = 2^2^^('^") 



exp 



00 



X'( 



2a 



:X(0o) + 20oX'(0o))e-2-s' + 



(61) 



(62) 



Thus, in the asymptotic region we have a behaviour exponentially close to a Randall- 
Sundrum behaviour, provided the parameter 7t is positive. This parameter is 



a 



K 



2M3 V 4M3y 



2 \ 1/3 



a/(e,0^/M3), 



(63) 



where we have introduced a function /(^, x) defined by 

/(e, x) = (^6 - C' - ^(2e - 1)) (1 - x^aY" 2Fi(1/2, 1/3, 3/2, ^x^A) 



+ r 



6 



X 



m - 1) 



(64) 



The coupling function is positive if ^0o < 4M'^. The positivity of k can always be satisfied 
with a large enough Brane-tension parameter a. Nevertheless, we can be more concrete 
by making a choice of 0(+O). We can take 0q = M^. Then, the positivity of the coupling 
function restricts the values of .^ to ^ < 4. The function /(^, 1) is negative for negative 
for all allowed values of ^. Thus, the warp factor will always be decreasing. The plot of 
this function of S, is shown in Figure 5. 

We can construct numerical solutions for the function A{y) and study the profile of 
the warp factor e^^"^' for different values of the parameters of the model. In this case, the 
brane tension determines the value of the derivative of A{y) at y = 0, so it fixes one of the 
two initial conditions needed for the numerical evaluation. The derivative of the brane 
tension with respect to the field, a', is proportional to a0o. Thus, changing it corresponds 
to a new value for 0(0). Note that, although we have restricted ourselves on 22-syKimetric 
solutions, asymmetric solutions are also possible. The profile of characteristic solutions 
is plotted in Figure 6. 

We see that, in general, the warp factor resembles the Randall-Sundrum decreasing 
exponential. Yet, for a range of values in the parameter space we get solutions which 
deviate slightly from this form. As the brane tension becomes smaller and ^ takes higher 
positive values, the warp factor exhibits a peak close to the Brane, before it starts 
decreasing again. This peak is amplified as we approach the value of ^ for which the 
coupling function /(0) tends to zero. 
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Figure 5: A graph of the function f(,^, 1) plotted for all values of ,^. THe function 
maintains a negative value. 
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Figure 6: The warp factor for different choices of A'[0), a, ^ and 0o in units of 

2M3. 

5 Smooth Spaces 



As we saw in the last example, the presence of the Brane was not essential to obtain a 
localized warp factor. In this section we shall consider solutions A{y) when the ''Brane" 
is the scalar field configuration itself. Such an example is well known in the ^ = case. 
Introducing a standard kink (0 = (pQ ta.nh.{ay)) into the .^ = equation of motion, we 
obtain in the Z2 symmetric case {A{0) = 0) 



3^(^) = (cosh(ai/))"^e-4t-°h2(aj/) 



(65) 



with 7 = 02/9M3. 

Smooth solutions of the Bulk equations of motion are also present for ^ 7^ 0. It is 
not difficult to see that the metric choice 



e^(y) = (cosh(ay))-^ , 
corresponds to the same scalar field solution (f){y) = 0o ta.nh.[ay) with 



7 = 2(r'-6), 0o = a-V(O) = (2M3)i/2, 



'6(1-60 
01-20 



(66) 



(67) 
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This solution exists for < .^ < 1/6 and only for the above specially chosen value of 0o- 

The curvature scalar of this space is R = 4a^7 f 1 — (1 + 57/4) l^- j . Note the asymptotic 

AdS value R{oo) = — 20a^(^~^ — 6)^. The scalar potential corresponding to this solution 
is a quartic function of the scalar field with tuned ^-dependent coefficients. Note that 
this solution is a particular case of (jHHJ). The above choice of 0o, together with the choice 
Co = 0, corresponds to 

2a 



Ci = , C2 = --^== v/^-i(l - 20(1 - 60 . (68) 

It is interesting that the same metric choice corresponds also to the solution 

(j){y) = 0(0) (coshiay))-' (69) 

with 02(0) = 12M^{^-^ - 6)/(3 - 160 defined in the same ^-range. 

Another intersting solution for the metric is defined by the choice 0o = 2a/^^^W^, 
for which ()58|) gives 

A = 2atanh(ay) i -2 + -(8 - C^) cosh-'^/^ {ay) 2^1(1/2, 1/3, 3/2, tanh^iay)) I . (70) 

Integrating, we obtain 

A{y) = -4 ln(cosh(ay)) + ^(8 - T') tanh^ay) Fppg ({1, 1, 7/6}, {3/2, 2}, tfxnh''iay)) , 



The warp factor e"^ is plotted in Figure 7 for ^ = 1/4 and for ^ = 1/9. 



(71) 




y 

0.5 1 1.5 2 2.5 

Figure 7: Warp factors for 0(y) = 0(0) (cosh(ay))~ , ^ < 1/2. 
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For larger ^ the behaviour does not change drastically. Note though that for values 
1^ > 1/2 the warp factor develops a maximum beyond the origin. In Figure 8 we plot the 
cases ^ = 1 and ^ = oo 




Figure 8: Warp factors for 0(y) = 0(0) (cosh(ay)) , ^ > 1/2. 



As already seen, the equation of motion for a kink-like scalar cannot be solved an- 
alytically for general boundary conditions. It is however possible to obtain numerical 
solutions. We expect to find a set of ^2-synimetric solutions for e^^^-* that reduce to the 
known ^ = solution mentioned in ^Sj. As an example, we consider numerical solutions 
of the warp factor equation, imposing the boundary values A{0) = 1, A'{0) = and tak- 
ing different values for 0o- The resulting warp factors for ^ = —2 and ^ = 0.8 are shown 
in Figure 9. We have taken (po = a = 1 in units of 2M^. For this choice of boundary 
values and units, ^ = 2 corresponds to the limiting value for which the function /(0) 
becomes zero at the origin, so higher values of ^ are forbidden. Notice the peak beyond 
the origin in the second plot. 





Figure 9: Smooth numerical solutions. 



6 Graviton Localization 

It would be interesting to check if gravity is localized in the geometries we calculated and 
especially those that deviate from the original Randall-Sundrum. Let us quickly provide 
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a general argument and then study some specific examples. We consider a perturbation 

SGmn = '^M'^iv^M^ (^' y)' (50 = , (72) 

for the gauge h^M = 0. Imposing transversality ( /i^ = d^h'^'^ = ), we obtain to first 
order, 

"^^ ^""^'^0' + A{y) + A' (y)] V = , (73) 



dy'^ / 

where d^ = rj^^d^d^ . Notice that this result is independent of the coupling function 
f{(p))- If we introduce a trial solution of the form of a product /i^,^ oc e^'^'^ipiii), we get a 
Schroedinger-like equation 

^' + A{y) + A\y)\ i;{y) = m^e-^^y^iy) , (74) 



dy 

where we have introduced the mass m^ = —p^. In order to study the spectrum of this 
equation, it is more convenient to transform it into a conventional Schroedinger equation. 
In order to eliminate the exponential, we may introduce the transformation 

^ = e-^^'^, ^{y) = e^/'^. (75) 

dy dz 



The resulting equation is 



with the potential 



d^ ,.,.>- 2 



2 + U{z) ) ^ = m'^lj , (76) 



dz 
^^, , 3rfM 9 (dAV 

Note that this equation can be put into the form 

d ?>dA\ ( d ?>dA\ - .- 

Tz-a-Tz)\Tz--aTz)'^ = '^'^- (^'^ 

This is supersymmetric Quantum Mechanics and the transformed graviton wavefunction 
(zero mode) corresponds to the supersymmetric ground state. This form also excludes the 
existence of tachyon modes. The zero mode is just ^0(2/) = Ne^^^^ and it is normalizable. 
We also have to know if there is a gap between the zero mode and the continuum of 
eigenstates. For this, we have to know the behaviour of the potential U{z). Although 
for most of the cases above the change of variable z = J dy e'^^^^^'^ is not analytically 
integrable, we may draw some conclusions with the help of the asymptotic behaviour 
dnH). Since, for y — > 00, we may have 

dye-^/'^ J dye^y/^g{y))-'/' (79) 

or z oc e^^/^. As a result, 

lim U{z) oc e^ ^ . (80) 

y-*oo 
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Therefore, the continuous spectrum starts from zero mass and there is no gap. 

Next, we may check the profile of the locahzation potential U{z) for various values 
of C,. For C = |, we have A{y) = —4:ln{cosh{ay)). In this case, the integration can 
be done analytically and the transformed coordinate is z = -^ sinh {2ay) + |. For ^ = 
0, the warp function becomes A(y) = — | (4 In (cosh (ay)) + tanh^ (ay)). We can only 
proceed numerically to perform the change of coordinates and calculate the potential. 
The resulting profiles are depicted in Figure 10. The localizing potential has the familiar 
volcano-like shape we also encounter in standard Randall-Sundrum. 



V(z) 





Figure 10: Localizing potentials for various values of ^. 

It turns out that the volcano-like profile of the potential is not maintained for all 
values of ^. For the boundary value and unit choices made for the smooth numerical 
solutions of section 5 depicted in Figure 9, we find that at the value ^ = (vl93 — 9)/16 ~ 
0.306 the global minimum at the origin y = changes into a local maximum. Thus, as we 
move towards higher ^'s, a central spike is developed. For ^ = 1/2 the potential becomes 
zero on the Brane, while as C, goes to infinity the potential at that point approaches 
unity. The corresponding graphs are shown in Figure 10. 



7 Conclusions 

In the present article we investigated the existence of solutions for a non-minimally 
coupled Bulk scalar field in a warped Brane-world framework. For a scalar field coupling 
to gravity of the form —^C^ip'^R, we derived a set of non-singular solutions for a wide range 
of the coupling parameter C,- We demonstrated the compatibility of the usual Randall- 
Sundrum warp factor with the presence of a non-trivial scalar field for a suitably chosen 
scalar potential. This was done for either a scalar field-dependent or independent Brane- 
tension. The profile of the scalar field solution in the field-independent Brane-tension case 
is that of a folded-kink ( tanh(a||/|) ). The scalar field acquires its minimum value on the 
brane, approaching a constant value as we move towards infinity in the y-direction. The 
conformal value of the coupling parameter C,c separates the above mentioned solutions 
from singular "solutions" . Thus, for the Randall-Sundrum warp factor, non-singular 
scalar field solutions exist only for C, > C,c and they are, in general, of the above folded- 
kink shape. For negative values of the coupling parameter ^ the scalar field "solutions" 
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found correspond to a scalar potential with negative powers of the scalar field and, 
therefore, singular. A field-dependent Brane-tension allows for a more diverse range of 
behaviour including scalar field solutions which exponentially decrease at infinity. 

Guided by the scalar field set of solutions with a folded kink type of profile, we 
investigated the existence of general warp factor solutions that are different from the 
exact Randall-Sundrum case but still localized. Thus, assuming a scalar field solution 
of the form tanh(a||/|), we derived corresponding warp factor solutions which we ana- 
lyzed semi-analytically and numerically. Our analytical treatment further showed that 
for a wide range of values in the parameter space of the model we get finite geometries 
which are well-behaved for large y, as long as the Brane tension we introduce is large 
enough. Furthermore, we considered smooth warp factor solutions in which the role of 
the Brane is played by the scalar field itself. In this setup we considered some special 
solutions and proceeded to study numerically general finite geometries. We concentrated 
on Z2-sjmm.ethc solutions, although asymmetric ones are also possible. We considered a 
class of solutions which asymptotically reduce to decreasing exponentials of the Randall- 
Sundrum type. These solutions exist for a coupling parameter C, within a range of values. 
For a subset of these localized solutions the warp factor is not a monotonous decreasing 
function but exhibits a second maximum close to but beyond the origin and subsequently 
decreases. We have also derived analytically special exact solutions existing for special 
choices of boundary values and for a range of the coupling parameter. For these solu- 
tions, the same warp factor corresponds to either a kink scalar solution or the solution 
= 0(0) (cosh(at/))^ . 

Finally we considered the localization of gravitons near the brane. Although, the 
Schroedinger-like equation for gravitational perturbations is the same as in the minimal 
case, the warp factor detailed profile depends on the coupling parameter and the details of 
the localized spectrum should depend on it. Of course, again the spectrum has no mass 
gap and does not contain any tachyonic modes. The form of the localizing (volcano) 
potential depends on the detailed profile of the warp factor. It was studied numerically 
in a number of cases but also analytically in special cases. For a particular choice of 
boundary scalar field value and the special coupling parameter value C = |, the localizing 
potential has the typical ^ = volcano profile. Nevertheless, for values of ^ larger than a 
certain value, the localizing potential develops a spike at the origin, that increases along 
with ^. For C, = ^ the spike reaches zero, while it tends to one for very large values of C,- 
This behaviour is currently being studied and will be the subject of a future publication. 
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